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Abstract 

We consider finite temperature properties of the Ising chain in a transverse 
field in the vicinity of its zero temperature, second order quantum phase tran- 
sition. New universal crossover functions for static and dynamic correlators 
of the "spin" operator are obtained. The static results follow from an early 
lattice computation of McCoy, and a method of analytic continuation in the 
space of coupling constants. The dynamic results are in the "renormalized 
classical" region and follow from a proposed mapping of the quantum dynam- 
ics to the Glauber dynamics of a classical Ising chain. 
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I. INTRODUCTION 



The finite temperature properties of systems in the vicinity of a zero temperature (T), 
second-order quantum phase transition have recently seen a great deal of attention |l] |J. 
In this paper we shall discuss properties of the quantum transition in an Ising model in a 
transverse field in spatial dimension d = 1. Apart from its intrinsic interest as a possible 
experimental system, the Ising model deserves scrutiny because it has a finite temperature 
phase diagram (see Fig |1] and its caption) which is extremely similar to that of the d = 2, 
0(3) quantum rotor model @||; the latter model is of relevance in the study of quantum 
Heisenberg antiferromagnets in d = 2. In Ref ||, a number of results were obtained for 
universal scaling functions of the d — 2, 0(N) rotor model by a 1/N expansion. In the 
present paper, we shall describe the corresponding scaling structure for the d = 1 quantum 
Ising model, and obtain exact results for crossover functions of some static and dynamic 
observables. 

An early study of the Ising model || obtained scaling results only at T = 0. Scaling 
results at finite T, but with transverse field set exactly at its critical value, can be obtained 
by conformal invariance (See e.g. Ref ||[7| and Section |11B| ). Other, previous finite T results 
include corrections to scaling at the critical transverse field ||, and results at T = oo which 
are beyond the scaling regime [(5, HI . 



In this paper we shall obtain new static and dynamic scaling results for the Ising model at 
finite T, with the transverse field away from the critical point; the ground state of the model 
is therefore "massive" . The results are for two-point correlators of the "spin" operator of the 
Ising model, and the static and dynamic results are obtained by rather different methods, 
(z) The static (or more precisely, equal time) results are obtained by using a lattice model 
computation in an early paper by McCoy pT| , and by a method of analytic continuation 
in the space of coupling constants. A separate scaling analysis of McCoy's results was 
also considered earlier in some limiting cases by Perk et. al. |§, who, unlike us, did not 
consider crossover functions. All but one of Perk et. a/.'s results agree with ours, and the 



disagreement is probably due to a miscalculation in the last section of Ref [fl l[]; this issue 
will be discussed in more detail in Section |11 A . 

(ii) The dynamic results are not rigorous, but are based on a physical argument which 
suggests that in the "renormalized-classical" region (Fig [I]), the long-time, long-distance 
spin correlations of the quantum Ising model should be described by an effective classical 
model. We propose that the appropriate effective model is the Glauber dynamics 0] of a 
classical Ising chain. This proposal then allows us to compute dynamic scaling functions in 
the renormalized-classical regime, up to a single unknown damping constant. 

The main technical innovation of this paper is the use of the method of analytic contin- 
uation in the space of coupling constants, which plays a central role in the computation of 
static correlations. This method was alluded to earlier by the author, Senthil and Shankar 
in Appendix C of a study [|J of the dilute Bose gas in d = 1; here we shall demonstrate 
its utility for the Ising model, and provide more details on the methodology for both the 
Ising model and the Bose gas. The method relies on the fact that one-dimensional quantum 
systems with short-range interactions do not have phase transitions at finite T. As a result, 
all observable properties should be smooth functions of the bare coupling constants at any 
non-zero T. There can be phase transitions at T = however, separating two phases with 
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very different physical properties — see the phase diagram and figure caption for the Ising 
model as function of coupling constant and T in Fig |l|. Raising T slightly above the two 
T = phases, gives finite T states whose properties are again very different. However, 
the observables in these states must be related to each other by an analytic continuation 
in the coupling constant. This relationship becomes somewhat less counterintuitive when 
one realizes that the two states have to connect through the intermediate "quantum-critical 
region" (see Fig |IJ) in which the observables only have a smooth, subdominant, dependence 
on the coupling constant ||. With the method of analytic continuation at hand, we can, 
from a knowledge of observables at low T above one T = phase, deduce the observable at 
low T above the critical point and the other T = phase. 

As a by-product of our method of analytic continuation, we shall obtain two apparently 
new integral representations of Glaisher's constant, Ac [0. This constant plays a central 
role in universal amplitude ratios of two-dimensional field theories whose correlators are 
described by Painleve differential equations |E,T3-rB|. We do not have analytical proofs for 



our representations of Ac, but will instead present numerical evidence which establish them 
beyond any reasonable doubt. One integral representation will appear in the analysis of the 
Ising model (Section |I A| and Appendix [A]), while the other will be deduced from the Bose 
gas results of Refs P|,|TT|] (Section |TTT| and Appendix |C]). 

Our main new results on the the Ising chain in a transverse field will be outlined in 
Section H, with computational details being provided in appendices [A] and [B]; the static 
results will be described in Section p A| , and dynamic results in Section p B| . In Section |II 
we shall review earlier results on the properties of the quantum transition in the dilute Bose 
gas in d < 2 |3J and display exact results for the universal crossover functions of some static 
observables in d = 1 as obtained in Refs |j,17,|]; some new results will also be established 
here, although the main purpose of this section is to observe that there are many strong 
formal similarities between these Bose gas results and the new Ising model results of this 
paper. Section [TV] will discuss some general issues, including the significance of the similarity 
between the Bose gas and the Ising model. This similarity will motivate a possible route 
to the solution of the remaining outstanding problem: the explicit computation of dynamic 
correlations in the quantum Ising model (with its unitary Hamiltonian time evolution), 
without recourse to a mapping to the phenomenological, dissipative, Glauber model. Such 
a direct solution, and its comparison to the phenomenological models, should increase our 
understanding of the relationship between the different approaches to dynamic phenomena 
in quantum statistical system. 



II. ISING CHAIN IN A TRANSVERSE FIELD 

This section will describe our new results on Ising chain in a transverse field. We will 
consider the Hamiltonian 

Hj = -jJ2 (9°x(i) + <Tz(T)<Tzii + 1)) (i) 
i 

where J > is an overall energy scale, g > is a dimensionless coupling constant, <r x (i), 
<j z (i) are Pauli matrices on a chain of sites i. This model has a T = quantum phase 
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transition || at g = g c = 1 from a state with long-range-order with (a z ) ^ (g < g c ), to a 
gapped quantum paramagnet (g > g c ). The critical exponents are z = 1 and u = 1. The 
crossover phase diagram at finite T is shown in Fig [I]. Notice the strong similarity between 
this and the phase diagram for the d = 2 quantum rotor model in both cases there 

is long range order only at T = for g < g c , and the physical interpretation of the finite 
temperature regions is very similar (see the figure caption of Fig []]). Further, at the smallest 
finite T above the ordered phase, the spin correlation length, £, diverges exponentially with 
T: 



kxt~(g c -gr/T. (2) 

In the d = 2 quantum rotor model this behavior of £ is due to the properties of the d = 2 
classical non-linear sigma model 0; in the d — 1 Ising model, the correlation length is 
of order the mean spacing between the thermally activated kink and anti-kink excitations, 
which have an energy ~ (g c — g) zv . 

We consider, now, the universal functions describing the crossovers in Fig |I] for Hi. 
First, it is necessary to define the appropriate renormalized couplings. It is known that in 
the vicinity of the critical point, the Ising model is described by a continuum field theory 
of free Majorana fermions |[L5|| . We choose the mass m and the velocity c of the Majorana 
fermions as our renormalized couplings (the energy of a fermion with wavevector k is there- 

1 /2 

fore (jn 2 c 4 + h 2 c 2 k 2 ^j ); notice that these renormalized couplings are defined at T = 0. The 
relationship between m and c, and the bare couplings of the lattice model is non-universal 
and, in general, difficult to determine; for the special case of Hj, which happens to be 
integrable, we have 

h 2 (g c -g) , 2Ja 

where a is the lattice spacing. The Majorana fermion mass m can have either sign, and 
moving m through tunes the system across the transition; we have chosen m > to 
correspond to the ordered side. Scaling would suggest that \m\c 2 ~ \g — g c \ zu , and the value 
of m for Hi agrees with zv — 1. The velocity, c, simply sets the relative distance and time 
scales and may be assumed to be constant across the transition. 

As we discussed in Section |, there can be no phase transitions in this Ising model at 
any non-zero T. We may therefore conclude that all properties of Hi are analytic functions 
of the bare coupling g at all non-zero T. Accidentally, because zv = 1, analyticity in g is 
equivalent to analyticity in m. This is a powerful principle which we shall use below. 

Universality appears in the continuum limit a — > at fixed m, c, and T (this requires 
J — > oo and g — ► g c ). An equivalent "condensed matter" perspective is that we are studying 
length scales much larger than a, energy scales much lower than J, but of order scales that 
can made out of combinations of m, c, and T . Following arguments in Ref 0, we may 
conclude that, in this limit 

(a z (x, t)a z (0, 0)) = Z{k B Tf/%i i-j-x, -j-t, _ , (4) 
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where Z is a non-universal normalization constant, and <&j is a fully universal crossover 
function; the above scaling form reproduces the T = results in the limit T \ 0. The non- 
universality of Z is linked with the anomalous dimension of the a z operator, indicated by 
the T 1//4 prefactor. This non-universality also means that a single normalization condition 
is necessary to set the overall scale of $/: this will be specified in the next subsection and 
leads to the value 

(5) 



Jl/4 



for the specific lattice Hamiltonian Hj. 

We emphasize that we are studying here correlators of the 'spin' operator a z . The 
expression for a z in terms of Majorana fermions involves an infinite string of fermion oper- 
ators ||15|| . In contrast, the 'energy' operator <j x is a bilinear in the Majorana fermions; its 
scaling properties are much simpler and easily computable JTj| . 



A. Statics 



This section will present our new results for the crossover functions of the two-point, 
equal time, correlators of a z . In particular, the long-distance limit of this correlator has the 
form 



( a z (\x\ -> 00,0)^(0,0)) ~ A! exp(-|z|/f/) 



(6) 



where Aj is an amplitude, and is the correlation length. The result (^) implies that these 
must obey the scaling forms 



-1 
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Aj = Z{k B T) l ' A gi 



mc 



k B T, 



(7) 



where fi(s) and gi(s) are completely universal functions of — oo < (s = mc 2 /k B T) < oo. 
One of our main new results is the explicit, closed-form, expression for these universal 
functions: 



ln#/(s) 



dy 



In coth 
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(8) 



where s can take any finite value between ±00, and 8 is the unit step function 

*>-{;;<?• m 
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Details of the computation are provided in Appendix |A]. We have set the overall scale of gi, 
and hence that of by the normalization condition 



lim 



gngj 

.sV4 



1. 



(10) 



A key property of the crossover functions fi(s) and gi(s) is that they are analytic for 
all — oo < s < oo. The apparent singularity at s = in the definition of fi(s) in (|8j) is 
deceptive: it is not difficult to show that all derivatives of fi{s) are actually finite at s = 
(see Appendix |A~|). The analyticity of gi{s) follows from its definition in (§), the analyticity 
of fi(s) and the fact that 



ds 



s=0 



(11) 



This last result can be obtained by an elementary quadrature, and ensures that there there 
is no singularity in the integration at y = in the definition of gi(s) in (H). Indeed, these 
analyticity properties were an important ingredient in the derivation of @ in Appendix [A|; 
there we take the continuum limit of a lattice model result of McCoy JllJ for s > 0, and 
obtain results for s < by analytic continuation in s. 

We show in Fig || a plot of the functions fi(s), gi(s). Some asymptotic limits can also 
be obtained by elementary means from (181) 



Ms) 




(12) 



The corresponding asymptotic limits for gi(s) are not so simple. In Appendix [A] we find 

( s 1 / 4 s — > +oo 
9l (s) = { & s = (13) 

s —> — oo 




The universal number C\ was computed numerically (Appendix [A]), and we found 



(14) 



where Aq is Glaisher's constant (given explicitly by lnAc = 1/12 — C'( — 1) where ((x) is 
the Riemann zeta function). We have numerically established that |e| < 10~ 20 , and there is 
essentially no doubt that in fact e = 0; if so, we have 



In (n^e 1 ^ 1 / 1 ^^) = J 



idy 

y 



' dfrjyY 

, d y . 



+ 



dy 



i y 



(15) 



as one of our new integral representations of Aq. We obtained our proposed value of C\ by 
using a conformal invariance argument to connect to a T = result of Pfeuty || (see Eqn 
(0) below). 
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Perk et. al. have also studied the continuum limit of McCoy's results [O for limiting 



regimes at low T, although they did not obtain crossover functions. The results in (|12l , |13l) 
are in agreement with theirs, except that their value for Aj at s = disagrees with our value 
for C\. The discrepancy at the critical point, m = 0, between their value for the amplitude 
Ai at low T and Pfeuty's || result at T = 0, led Perk et. al. |J to suggest that the crossover 
function connecting to the T = result (Eqn (|lj]) below) must break down. Our results 
show that there is nothing wrong with the crossover function in fllEp (indeed, it is demanded 
by conformal invariance), but it is the value of Aj at s = in Ref || which is incorrect. 
The error appears to be traceable to an earlier error in the Section 6 of McCoy's paper BTT 



which computed correlations at m = and finite T. In the present paper we have only used 
McCoy's results for m > from his Section 3, and have verified his computations in this 
case; for completeness, Appendix [A| contains a sketch of the required intermediate steps. 



B. Dynamics 

We begin by reviewing and remarking on some earlier results for dynamic correlations. 
(i) The scaling function can be obtained exactly at m = by using the conformal 
invariance of the T = critical theory. One finds |TS|,P,|r9|^JSO in imaginary time: 



2- 1/8 ^(0) 
(cosh(27rx) — cos(27tt)) j 



^(x, -ir, 0) = ; _ _ x (16) 



where we are using x = (kBT/hc)x and r = (ksT /K)t. Note that 

^(1x1-0,0,0) = ^,, (17) 

which agrees with the T = result of Ref , while 

$ / (N^oo,0,0) = ( ?/ (0)e-^ l/4 , (18) 



which agrees with (p|,pT3|) ) ; this is the agreement between the low T and T = results claimed 
earlier. We can also Fourier transform ([TBI) an d analytically continue to real frequencies to 
obtain the dynamic susceptibility x{k, w), which is the Fourier transform of the retarded a z , 
a z commutator. There are a number of interesting subtleties that appear in this process — see 
Refs here we display the final result: 

(1 .uJ + k\ (1 .uJ-k\ 

, (k ::)- Z9l{m7/8) hC ' ^ > 1 I MCA 

X[K ' UJ) 47rr(l/8) (k B T)V* (15 W+k\ (15 u-kY 1 J 

We have used k = hck/ksT and u; = hu/kBT; recall also that this result is valid only at the 
critical coupling m = 0, g — g c . A plot of the spectral density Imx(fc,w)/w as a function of 
lo is shown in Fig |3| for a number of values of k. Notice that there is relaxational behavior 
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for small k (e.g. k = 0.5) with a peak in the spectral density at uj = 0; at large k the peak 
moves to u ~ k indicating a crossover to the critical excitations of the T = ground state. 
(ii) Perk et.al. || obtained a set of non-linear, partial differential equations satisfied by an 
observable related to $/. However, T does not appear as a parameter in any of the equations, 



which are therefore identical to those obtained at T = (see Ref [15 and references therein). 

All effects of temperature appear in the initial conditions, which presumably have to be 

specified on some space-like surface. At the critical transverse field, Perk et. al. || used the 

differential equations to conjecture the form (|16|), apart from the overall scale (recall the 

discussion in Section |11 A| on problems in Ref || with the overall scale). It remains to be 

seen whether this approach can yield results on away from criticality. 

(Hi) A number of studies have examined finite T, time-dependent, on-site (x = 0) correlations 

of the lattice model at the critical coupling g — g c ||. The scaling limit of all of these results 

is contained already in (|TBJ); instead, they contain a great deal of additional information on 

the non- universal corrections to scaling in the nearest-neighbor model Hj. 

(iv) We also recall results of Perk et.al. PJTO]] for time-dependent correlation functions at 



T = oo for the lattice model. These results are dominated entirely by non-universal lattice 
effects, and contain no overlap with our universal scaling form @. 

The major gap in our knowledge of time-dependent correlations is the form of the dy- 



namic scaling function away from the critical-coupling (m ^ 0). In Section [IV|, we will sketch 
a strategy by which these may be obtained. In the remainder of this section we will use 
some phenomenological arguments to obtain some new results on the form of the dynamic 
scaling functions in the renormalized-classical region of Fig. [[]. 

Notice that in the renormalized-classical (RC) region (mc 2 3> k R T), the correlation 
length becomes (from (|8[), flT2| ) and also Ref ||) exponentially large: 




^(asbfl exp l^l (20) 

In particular, is significantly larger than the length scale ~ h/mc at which critical 
fluctuations are quenched, and the ordering of the T = ground state appears. Using 
arguments similar to those for the d = 2 rotor model [fj], we may conclude that at length 
scales ^> h/mc, but of order £rc, an effective classical model of the system should be 
adequate. We propose here that the appropriate classical model is the Glauber dynamics 



of the d = 1 classical Ising chain | I2| . Formally, the above arguments mean || that in the 
RC region, the scaling function, $7, of 3 arguments in (f|) collapses into a reduced scaling 
function, (fine, of 2 arguments: 



x Tit \ Ti 
(a z (x, t)a z (0, 0)) = N'^ R c ( — , rnc 2 > k B T, \x\ > — , \t\ > -— , (21) 

Arc m£ RC J mc mc 2 



where No = Z l / 2 [mc 2 ) 1 ^ = (a z ) T=0 is the ground state magnetization, and 4>rc is precisely 
the dynamic scaling function of the low T limit of the Glauber model |12[ (a very similar 
collapse occurs in the Luttinger liquid region of the dilute Bose gas; see Sec |TTT| and Ref. J3J). 
Notice that the time t scales as £ RC — this is because the Glauber model has dynamic ex- 
ponent z = 2. The dimensionless constant 7 determines the relaxation rate constant of the 
Glauber model; on dimensional grounds we expect 
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7 = «(^Y (22) 



mc 2 



where 7?. is a universal number, and p an unknown exponent. The functional form of the uni- 



versal function <Prc can be determined by taking the scaling limit of Glauber's solution [12 
at low T; we obtained (see Appendix [B]): 

i; \ „— f I I I-Ti ~ \ 

(23) 




where erfc is the complementary error function. Notice that (f)Rc(f, 0) = e ' r '; combined 
with (|2T|) , this result agrees with the mc 2 3> &bT limit of (|4])-(1l3[). 

III. DILUTE BOSE GAS 

In this section we will review earlier results, and obtain some new results, on the finite T 
crossovers in the dilute Bose gas in d < 2. There are two main reasons for doing this here: 
(i) We will show that the results for the crossover functions have strong formal similarities to 



the corresponding quantities for the Ising model. In Section [TV], this similarity will motivate 
proposals for computing dynamic correlations directly in the quantum Ising model. 
(m) There are two distinct forms for the crossover functions in the literature — one obtain 
by Korepin et. al. H0], and other in Ref. [|J by the method of analytic continuation. In 
this section, and Appendix 0, we will establish that these two forms are in fact equal, and 
in the process obtain a second new integral representation of Glaisher's constant. 

We begin by reviewing known results for the Bose gas. Consider a dilute gas of non- 
relativistic, spinless, bosons of mass m#, at a chemical potential fi, in d spatial dimensions. 
The partition function of the system is given by 



r ( 1 r n/k B T 

Z = yD*exp --^ drC(r) 



C(t)= I ' d d x M!\ x ,t)?^P^- - — ¥(x,t)V 2 ^(x,t) -fi\^(x,r)\ 2 
J [ or 2m b 

+ d d xd d x'\V(x,T)\ 2 v(x - x')\V(x',t)\ 2 , (24) 

where *&(x, r) is a canonical Bose field over spacetime x, r. The bosons experience a re- 
pulsive, short-range, but otherwise arbitrary, mutual interaction, v(x). At T = 0, there is 
a continuous quantum phase transition as a function of /i at // = 0, when the density of 
particles in the ground state changes from zero to a finite value. This transition has dynamic 
exponent z = 2 and correlation length exponent v — 1/2 [p0[| . The finite T crossovers near 
this quantum critical point are especially universal for d < 2 (in the sense that everything, 
including the scale factors, is independent of the details of the interaction v(x)), and are 
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summarized in the phase diagram in Fig f| ||. In the vicinity of the critical point, the 
two-point correlator of ty(x,t), where now t is real time, obeys 0] 



V(x,t)&(0,0) 



2m B k B T 

~~h 2 



d/2 



^ / y/2m B k B T k B T ^ /i 



(25) 



with $b a completely universal scaling function, independent of the boson interaction po- 
tential. Note that, unlike the Ising model, there are no anomalous dimensions and it was 
not necessary to define any renormalized couplings. Only bare coupling constants appear in 
the scaling form: this "no-scale-factor universality" is a consequence of the fluctuationless 
ground state for /i < 0, and the ultraviolet finiteness of the critical quantum field theory [|J. 
These last properties are rather special and, in particular, are not present in the Ising model 
considered below. In d = 1, closed-form expressions for crossover functions, corresponding 
to some limiting cases of (|2~5"D , have been obtained ||[17],[|]. The long-distance limit of the 



equal time correlator is argued to have the leading behavior 



#(|a;| 00,0)^(0,0)) ~ A B exp{-\x\/£ 



IB 



(26) 



From (|25|) we see that the amplitude A B and the correlation length £ B must satisfy the 
scaling forms 



! _ y/2m B k B T ( /i 

<5B — Z JB 



h 



k B T 



. \J1m B k B T ( fi 
a b = g B 



h \k B T , 

The crossover functions, f B , g B , in the form obtained in Ref M are 



(27) 



f B (s) 



00 dy , , \y 2 — s r—, nl 

— lncoth^ l - + J\s\9(-s) 

on 2 v 



ln# B (s) = -2^ dy y 



dy 



lnC 2 



(28) 



where C2 = tt ' 2 e ' 2 2 ' 3 Aq , and —00 < (s = jj,/k B T) < 00. As in the Ising case, both 
crossover functions are analytic functions of s, for all — 00 < s < 00. In particular, and 
despite appearances, the result for f B (s) in (pg|) is analytic at s — 0. Refs |6|,|17|| give 
different expressions for the crossover functions in the region s < 0, whose equivalence to 
( p8|) has not so far been explicitly established — we provide the missing proof in Appendix y. 
The proof requires a new identity involving Glaisher's constant very similar to, but distinct 
from, that obtained in Section [11 A| : 



In (tt 1 /^ 1 ^ 2 / 3 ^ 6 ) =2 f dy 



' df B (yY 

, d y . 



+ 



Mj 



+ 2 



-1 



dy 



' df B (yY 

- dy , 



(29) 



Again, we do not have an analytical proof of this identity, but have obtained overwhelming 
numerical evidence for its validity (Appendix [C]). 
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The striking formal similarity between the results ( p8|) for the quantum transition in the 
dilute Bose gas, and the results @ for Ising model should be readily apparent. Transforming 
to dimensionful quantities, the main difference between the results for the correlation length 
is that the non-relativistic dispersion U 2 k 2 / (2m#) — [i of the z = 2 Bose gas has been replaced 

by the relativistic (m 2 c A + h 2 c 2 k 2 ^j in the z — 1 Ising model; otherwise the functional 
forms are essentially identical. The expression for gs in terms of Jb is also similar to that 
for gj in terms of //; the latter expression contains an additional integral / dy/4y, which is 
surely related to the presence of anomalous dimensions in the Ising model. Note, however, 
that the asymptotic behaviors of //, gi are quite different from those of gs- 

To emphasize the significance of the above similarity, we note that it is present despite 
some important differences between the Bose gas and the quantum Ising model. Both models 
are 'solved' by mapping to models of free fermions — however, the Bose field ^ b only has 
non-zero matrix elements between states whose fermion number differ by unity; in contrast, 
the Ising spin a z has non-zero matrix elements between states with arbitrary numbers of 
fermions. Second, the \i > ground state of the Bose gas has gapless excitations, while the 
Ising model has a gap for all m ^ 0. 



IV. DISCUSSION 

In this paper, we have obtained universal crossover functions (Eqn. @) of static ob- 
servables of the finite T quantum Ising chain. The crossover functions turned out to be 
remarkably similar in form to those obtained earlier for the dilute Bose gas P,|TT|,p| . This 



similarity suggests that the sophisticated mathematical structure described in Ref to 
describe correlators of the Bose gas, should also apply to the continuum limit of the Ising 
model in a transverse field. Application of these methods to the Ising model should yield 
more information on the scaling function $/; in particular we should be able to study un- 
equal time correlators away from the critical point. Traditional analyses of the Ising model 
have been carried out on a lattice |5| JTT|j25|| , but the simplifications suggested here appear 
only in the continuum limit. It would clearly be preferable to use, instead, an analysis 
which is valid in the continuum limit from the outset. Fortunately, just such an analysis 
has appeared in some recent work [21, IB]; these papers use non-trivial results on the form 



factors of the continuum Ising model |22j to deduce properties of the correlation functions. 
The existing results are at T = 0, but the author has some preliminary results on extending 
the determinant representation of Ref [16j to non-zero T ||23|| . Application of the methods 



of Ref H to this determinant representation is the next step, but has not been carried out 
yet. 



Such an analysis should be able to test the conjectured RC scaling forms of Section [II B 
by a direct solution of the quantum problem. If achieved, such a result would yield the exact 
value of the dimensionless 'rate constant' 7. Further, it would make a connection between 
two rather different approaches to dynamic phenomena — the quantum dynamics couples the 
Ising model to a heat bath only to set the initial density matrix exp(— Ti/ksT), and the 
subsequent evolution is completely unitary; in contrast, the classical Glauber dynamics has 
the spins in constant contact with a heat bath. Notice that the classical Glauber dynamics 
has purely relaxational behavior for the a z spins; such a behavior might be questioned on 
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the grounds that Hi is integrable and possesses an infinite number of conservation laws. 
However, the conservation laws are expressed simply using the fermionic degrees of free- 
dom, and are highly non-local in terms of the a z variables; as a result, we believe that the 
integrability will not disrupt the expected mapping to the relaxational Glauber dynamics; 
for a related discussion on the role of integrability, see Ref |f2^}] . In any event, it is clear 
that a direct quantum study of the dynamics away from criticality will surely enhance our 
understanding of dissipative dynamics in macroscopic quantum systems. 
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APPENDIX A: DERIVATION OF CROSSOVER FUNCTIONS OF THE ISING 

MODEL 



It appears worthwhile to sketch a derivation from first principles using a consistent 
notation. 

First, following Lieb, Schultz and Mattis ||, convert Hi into a free fermion Hamiltonian 
by the Jordan- Wigner transformation. Then, evaluate the equal-time spin correlator in 
terms of the free-fermion correlators. This yields an expression for the correlator in terms 
of a Toeplitz determinant p|,p5|: 



(vz{i)&z{i + n)) 



Go G-i 
G 1 



G- n 



+ 1 



G 



n-l 



Go G_i 
G\ Go 



(Al) 



where 



with 



G(4>) 



G„ 



1-ge* \ 1/2 



2tt 



27T 



1 - ger 



tanh 



J 



knT 



((1 - ge «)(l - ge-«)) 



i<M\l/2 



(A2) 



(A3) 



We can now take the large n limit of (|A1|) by Szego's lemma [13[] provided g < g c = 1. For 
g > g c a more complicated analysis is necessary to evaluate ( |A1| ) directly ||11|| ; we shall not 
need this here as we shall obtain results for g > g c by our method of analytic continuation. 
By Szego's lemma we obtain 
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fim (a z (i)a z (i + n)) ~ e nA ° exp [ ^ pA p A_ p ] , g < 1 



(A4) 



where 



lnG(0) = £ V* 



(A5) 



The correlation length is clearly fa = —a/X , and the above approach gives us the value 
of fa for g < g c . From (|A3| - |A5|) we can deduce that 



In coth 



J 



k B T 



1 + g 2 — 2g cos ( 



,1/2 



, 9 < 1 



(A6) 



Now change variables of integration to = 4>/a, re-express in terms of renormalized couplings 
by writing J = Tic/2a and g — 1 — mca/h, and take the scaling limit a — >• 0. This yields 



00 d& 
2^ 



In coth 



m 



2 C 4 + /j2 c 2 fc 2W2 



m > 



(A7) 



Let the scaling function for the inverse correlation length, in (0), //(s) = /^(s) for s > 0, 
where recall that s = mc 2 /ksT. Then clearly 



1 

7T 70 



rfy In coth 



(y 



2 , s 2Xl/2 



(A8) 



Notice that //"(s) is an even function of s. Later, it will be useful to have another form for 
ff"(s): change variables of integration y — > (y 2 — s 2 ) 1 ^ 2 , and integrate by parts to obtain 



7r J\s\ smh-y 



| sinh y 



(A9) 



We now need to get the scaling function of the inverse correlation length, fi(s) for s < 0. 
We could get this by returning to an analysis of (|A"T|), but as noted earlier, Szego's lemma 
does not directly apply and the analysis is not straightforward ]TI |. Instead, we use the 
method of analytic continuation to obtain the answer quite rapidly Unlike the needed fi(s), 
the function //"(s) is not analytic at s — 0. Let us examine the nature of the singularity in 
/f(s) near s = 0. Insert 1 = ((y 2 + s 2 )/(y 2 + l)) 1/2 ((y 2 + l)/(y 2 + s 2 )) l/2 in the logarithm 
in (Kgh to set 



1 r°° 
f>(s) = - / dyln 

7T JO 



v 2 + l 



coth 



(y 2 + s 



2U/2 



+ ^ / */ ln ( y 2 r 1 1 



(A10) 



The first integral has an integrand which is a smooth function of s 2 for all values of y 
(xcothx is regular at x = and has a series expansion with only even powers of x); so 
for small s it can only yield a power series containing even, non- negative powers of s. The 
second integral can be performed exactly and we obtain 
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f^(s — > 0) = — |s|/2 + terms with even, non-negative powers of s (All) 

The — \s\/2 term is the only non-analyticity at s — 0. It is then evident that defining 

fj(s) = f>(a) + \s\9(-s) (A12) 

gives the unique analytic function which equals ff (s) for s > 0; this was the result quoted 
in (H|) for the crossover function of the correlation length, and is valid for all values of s. 

We turn next to the scaling function, gi(s), of the amplitude. We need to evaluate the 
summation in 



techniques for doing this were developed in Ref |TT| and we shall not 
reproduce them here. The result is contained in Eqn. (3.27) [^6| of Ref PI , and taking 
its scaling limit in a manner similar to that discussed above for the correlation length, we 
obtain 



kigi(s) 



2tt 2 
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sinh y J s sinh y 
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,2)1/2 



,2)1/2 
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,2)1/2 
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r,2 



,2)1/2 



+ - Ins 
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s > 0. (A13) 



We now have to analytically continue this result to obtain the scaling function for all values 
of s. We first take its derivative and obtain 



9i[s) 



ds 



7T Z 



dy 



sinh y 



dy s 1 



Remarkably, the integrals over y and y have factorized. Now notice by comparing with 
that 



1 dgijs) _ 
9i(s) ds 

From ( |A13j ) it is evident that In gi(s 
may integrate (|A15 ) to obtain 



ds 



(A15) 

(1/4) Ins; using this boundary condition we 



9i{s) 



exp 



dy f df>( y y 

y v d v , 



'dy 
Ay 



s > 



(A16) 



Because df^ / dy\ y \ = —1/2, the above integral is a smooth function of s as s \ 0. The 
result (H) for (?/(s), for all s, now follows simply by analytically continuing ff(y) to fi{y). 

Finally, we consider the asymptotic limits in (|T3|). Those for fi(s) can be obtained by 
elementary methods. The limit gi(s — > +oo) is also evident from (||). Consider next gj(s — > 
— oo). Use (|5p, and express fi(s) in terms of /^(s); after some elementary manipulations 
we obtain 



ln5f/(s 



-oo 



In I s I — 



l dy 
o y 



r 2 dfm + ; 

v dy 



dy dff{y) 
i y dy 



(A17) 



After using ( |A8| ) , and interchanging orders of integration, all integrals can be found in 
standard tables. The integrals in (|A17| ) actually sum to 0, and we obtain the result in (|13|): 
lng/(s — > — oo) = — (3/4)ln|s|. 
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Finally, consider the value of gi(0). We have been unable to evaluate the integrals in 
(§) analytically. Instead, we carried out a high precision numerical integration, using the 
arbitrary precision numerics built into the Mathematica package. This gave us 

In 0/(0) = -0.152318694592340635 (A18) 

Compare this with the numerical value of the expression obtained by combining the T = 
results of Pfeuty with conformal invariance: 

In (7r 1/4 e 1/4 2 1/12 A G 3 ) = -0.152318694592340634987799466 .... (A19) 



There is essentially no doubt that the expression in ( |A19 ) is the exact value of In 0/(0), and 



this gives us one the claimed representations of Glaisher's constant. 

APPENDIX B: SCALING LIMIT OF GLAUBER'S SOLUTION 

I thank S. Majumdar for showing me how to take the scaling limit discussed in this 
Appendix. 

Glauber [|nj considered the classical Ising model (Eqn ([TJ) at g = 0) coupled to a heat 
bath at a temperature T. He derived a differential equation for the correlator C(i,t) = 
Mi,t)<r z (p,0)): 

= ~ C (h t) + l tanh t^F - M) + C(i + 1, t)} (Bl) 
a at 2 kb± 

where a is a relaxation rate constant. At low T, the classical Ising model has a correlation 
length £ = (a/2)e 2J / fesT (a is the lattice spacing) and the equal time correlator obeys C(x = 
ia, 0) = e~' x '^. As £ ^> a, we may take the continuum limit of ( [B~T| ) and obtain the partial 
differential equation 

1 dC(x,t) ld 2 C(x,t) C(x,t) 



aa 2 dt 2 dx 2 2£ 2 

This equation is to be solved subject to the initial condition C(x, 0) = e~' x '^. This is easily 
done by a spatial Fourier transform and we obtain the result in (^1|,^), with the unknown 
constant 7 proportional to aa 2 . 

APPENDIX C: COMPUTATIONS FOR THE BOSE GAS 

In this appendix we will establish the equivalence of our results for the crossover functions 
of the dilute Bose gas, /b(s) and 0s (s) in (|2"8"D, to those for s < in Refs |T7[ . 

As in the case of the Ising model in Appendix it is useful to define the function fs(s), 

f>(s)= r^lncoth^f^, (CI) 
Jo n 2 
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such that /b(s) = + y/^s9(— s). The proof that is analytic at s = is very 

similar to that presented in Appendix [A] for fi{s), and will be omitted. 

Let us call gsK, the result for the amplitude gs{s) obtained in Ref for s < 0. 



Then |7| 



9bk{s) 



exp 



/^lncoth^ + 2 



tt Jo y 



dy 



s < 



(C2) 



Clearly we have lim^^oo \J — ^sgBK{s) = 1. For lim s ^_ 00 ^—Asgsis) to equal unity, we 
must have from (1281) that the universal number 



dy 




-l V 



(C3) 



(C4) 



is equal to 

In (7r 1/2 e 1/2 2 2/3 A G 6 ) = 0.04193620109529 . . . 

We evaluated (|C3| ) numerically and found 

C 3 = 0.0419362011 (C5) 

which is in complete agreement with (|U4]) . The equations ( |03| , |C"4] ) therefore contain one of 
the claimed integral representations of Glaisher's constant, Ac- The integrand in ( |C3| ) only 
has a power-law convergence at +oo, and hence the numerical integration is not as accurate 
as that leading to (|A18|) , which converged exponentially fast. This difference between the 
two representations is directly linked to the fact that the Ising model has a gap at T = 
for all m ^ 0, while the Bose gas has gap less excitations for > 0. 

Let us now assume that linis^^oo ^—Asgsis) = 1. Then proving the equality between 
(7s(s) and (?bk(s) for all s < can be shown (after substituting /b(s) = fs( s ) + 9(—s)V—s 
in Q2"8|) and performing some elementary manipulations) to be equivalent to establishing that 



dy df> 



B 



— S 



t=y dy 



71 



dy 



In coth 



y 2 - s 



s < 0. 



o y 



(G6) 



As both sides of the equation vanish as s — ► — oo, it is sufficient to establish the equality 
of their derivatives with respect to s. Taking the derivative, and then converting all of the 
In coth terms to 1/ sinh by integrating by parts, one can show after some straightforward, 



but tedious, manipulations that (|Cq ) is indeed true. 
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FIG. 1. Crossover phase diagram for the d = 1 Ising chain Jij (Eqn. (p])). The quantum phase 
transition at T = 0, g = g c has exponents z = 1, v = 1. There is true long-range order at T = 
for g < g c , as indicated by the thick line. The parameter m ~ g c — g, so m = at the critical point 
and m > on the ordered side. An essentially identical phase diagram applies to the d = 2, 0(3) 
quantum rotor model P||, with the difference that the exponents z = 1, u ~ 0.7; the crossover 
phase boundaries, which occur at \g — g c \ zv ~ T, would then not be linear. All regions display 
critical fluctuations at short scales. In the renormalized-classical and quantum-disordered regions, 
these cross-over to behavior appropriate to the non-critical ground state before thermal effects 
have become apparent. In contrast, in the quantum-critical region, the critical fluctuations are 
quenched in a universal manner by thermal effects, and characteristics of the non-critical ground 
state do not appear at any scale Q. 
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FIG. 2. The universal crossover functions fj and gi of the Ising model as a function of 
= mc 2 /ksT 
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FIG. 3. The spectral density of the dynamic susceptibility, Imx(k, lu)/u, of the transverse 
field Ising model in d = 1 at its critical coupling g = g c (m = 0) but with T finite. The scale of the 
vertical axis is determined by the non- universal constant Z in (|l^), and has been chosen arbitrarily. 
We have k = hck/ksT, UJ = huj/kBT; the horizontal axis is therefore measured in dimensionless 
units and its scale is meaningful. 
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FIG. 4. Crossover phase diagram for a dilute Bose gas in d = 1 from Ref @ . The T = quantum 
phase transition occurs at chemical potential \x = and has exponents z = 2, v = 1/2. The 
temperature T is the most important energy scale in the quantum-critical region. The Luttinger 
liquid has power-law correlations at T = (indicated by the thick line), while there is exponential 
decay at any finite T. The same phase diagram applies for 1 < d < 2 except that there is true 
long-range order at T = on the thick line, and the nomenclature "incipient superfluid" is more 
appropriate than "Luttinger liquid" . 
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